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Introduction 

Me  will  consider  some  tuo   diraenslonal  flows  of  an 
Incompressible  heterogeneous  ideal  fluid  of  finite  depth  with  a 
free  surface.   Let  the  coordinate  axes  be  chosen  so  that  when  the 
fluid  is  at  rest  it  occupies  the  strip  -h<y<0,  -  ok  x<  oo,  of  tjfie 
x-y  planej  the  line  y  =  0  corresponds  to  the  undisturbed  free 
surface,  and  y  =  -h  to  the  rigid  bottom.   The  density,  p,  of  the 
fluid  when  at  rest  is  assumed  to  depend  on  y  only,  and  in  order  for 
the  equilibrium  to  be  stable,  p  =  p(y)  must  be  a  monotonic  non- 
increasing  function  of  y.   The  fluid  is  then  said  to  be  stably 
stratified. 

The  study  of  heterogeneous  fluids  was  begun  by  Stokes  (I8I4.7) 
who  considered  infinitesimal  gravity  vjaves  in  a  fluid  consisting  of 
two  layers  of  different  density.   This  was  later  generalized  to  any 

niomber  of  layers  by  VJebb  and  by  Greenhill  (see  [2  J  §a  231,  233 )« 

For  a  fluid  of  n  layers  w'.th  a  free  surface,  this  problem  leads 

to  the  following  conclusions.   There  are  n  "critical"  velocities, 

*-•<  U^  1  <  U  o  ...<  U  <  CO  :foJ^  snj   velocity  between  0  and  U  , 
n—x    Ti—d'  o  n-x 

there  are  n  possible  waves;  betx;een  U  t  and  U  _2  there  are  n-1 
possible  waves;  and  so  on,  vuitil  for  U  >  U  no  waves  are  possible. 
The  shortest  wave  for  any  given  speed  is  similar  to  the  surface 
wave  in  a  homogeneous  medium.   All  the  other  vjaves  are  internal 
waves,  having  one  or  more  nodes  below  the  free  surface.   As  the 
velocity  approaches  one  of  the  critical  speeds  from  beloxij,  the 
length  of  the  longest  wave  Increases  monotonically  to  00 .   Some 
interest  in  the  problem  derives  from  the  fact  that  under 
appropriate  circumstances  internal  gravity  waves  can  be  very  large 
and  their  effect  quite  noticeable,  as,  for  example,  in  the  case 
of  "dead  water"  resistance  (see  [2],  §§  231,  21^9). 
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In  contrast  to  the  multilayered  case,  the  problem  of  gravity 
waves  for  a  fluid  with  a  continuous  density  distribution  has  not 
been  treated  extensively.   Although  the  equations  were  derived  by 
Love  in  I89I  [1;]  only  a  few  special  cases  have  been  investigated 
(see  [2],  I  235,  [3],    [6],  and  references  quoted  there).   Current 
interest  in  the  problem  is  in  part  due  to  the  possible  bearing  on 
atmospheric  motions. 

This  report  is  concerned  with  gravity  waves  of  small  amplitude 
in  a  fluid  with  a  continuous  density  distribution.   For  the  most 
part  it  will  be  assumed  that  p(y)  has  a  piecewise  continuous 
derivative,  p'(y),  and  that  p'(y)  <  0,  but  is  otherwise  arbitrary. 
Section  I  deals  with  some  of  the  general  properties  of  progressive 
waves  in  the  medium.   This  leads  to  a  boundary  value  problem  for 
an  ordinary  differential  equation  with  two  parameters,  one  of  which 
also  enters  into  the  boundary  condition  at  the  free  surface.  As 
in  the  multilayer  case,  one  concludes  that  there  are  a  finite 
number  of  possible  waves  for  any  given  velocity.   However,  the 
number  is  not  limited,  since  there  are  an  infinite  number  of 
critical  speeds.   Again,  the  shortest  wave  behaves  in  a  different- 

way  from  all  the  other  x^Iaves,   It  is  proved  that  there  cannot 

1/2 

exist  waves  with  a  speed  greater  than  (gh)  '  ,  the  critical  speed 

for  a  homogeneous  medium.   An  example  is  given  of  a  fluid  of 
infinite  depth  which  also  exhibits  an  infinite  number  of  critical 
speeds. 


In  Section  II  we  examine  the  steady  state  problem  for  a 
surface  pressure  moving  with  a  constant  speed.   An  integral 
representation  for  the  perturbed  stream  function  is  derived,  and 
the  behaviour  of  the  flow  far  away  from  the  pressure  disturbance 
is  investigated.   A  formula  for  the  wave  resistance  is  given,  and 
from  it  the  conditions  for  the  vanishing  of  the  resistance  are 
deduced. 

The  uniqueness  question  for  the  steady  problem  suggests  the 
study  of  the  time  dependent  solution  due  to  a  moving  surface 
pressure  (see  [?],  pp.  210-218),   The  solution  for  general  p(y) 
is  given  in  Section  III.   However,  the  proof  that  this  solution 
approaches  the  steady  solution  is  done  only  for  the  special  case 
where  p  =  e~P^  . 
Section  I.  Waves  with  a  free  surface 

A.   The  equations  for  the  motion  of  a  two  dimensional  incompressible 
fluid  are: 

(1)  p(u^  +  UU^  +  VUy)  +  p^  =  0 

(2)  p(v^  +  uv^   +  vvy)  +  Py  +  gp  =  0 

(3)  P^  +  (pu)^  +  (pv)   =  0 

(^^      It  =  Pt  -^  ^Px  ^  ^Py  =  ^ 

Here  u  and  v  are  the  velocity  components  in  the  x  and  y  directions, 
p  is  the  pressure,  and  p  -  the  density.   Equations  (1)  and  (2)  are 
the  Euler  equations  of  motion,  (3)  is  the  equation  of  m.ass 
conservation,  and  (L|.)  is  the  incorapresslbility  condition.   In  view 


^  d) 


of  ik) f    (3)  can  be  replaced  by  the  more  convenient  relation 

(5)  ^x  "^  ^y  =  °  • 

The  kinematic  free  surface  condition  is 

^^^      ^=       %■*"  "^x  =  ^        ^*      y=f^(x,t). 

If  a  surface  pressure,  P(x,t),  is  applied,  the  dynamic  condition  is 

(7)  p(x,>7(x,t),t)  =  P(x,t)  . 
At  the  bottom  vje  have 

(8)  V  =  0     at    y  =  -h  . 

If  P  =  0  and  the  fluid  is  in  a  state  of  rest,  the  solution  is 
u  =  V  =  0,     "^(=0,      p  =  g  V  p(s)ds 

y 

where  p  =  p(y)  is  continuous  and  monotonic  decreasing  with  y.   The 
equations  for  a  small  disturbance  from  the  state  of  rest,  brought 
about  by  the  application  of  a  small  surface  pressure,  can  be 
obtained  formally  by  the  standard  pex^turbatlon  procedure.   Since 
equation  (5)  for  the  perturbed  quantities  does  not  change  in  form, 
we  can  introduce  a  stream  function,  ^,  for  the  perturbed  velocities. 
We  will  designate  these  velocities  by  u  and  v  again,  since  no  other 
velocities  will  be  referred  to  from  this  point  on.   Thus,  take 
\j((x,y,t)  such  that 

We  can  eliminate  the  perturbed  pressure  and  density  from  the 
equations,  and  finally  obtain 


f^rfrt 


:\    U   ^-c!' 


(9)  [(P^x)x   +    (P%)yUt    -    SP'^XX   =    0 

where    p  =    p(y)    is   the   originally  prescribed  density  distribution. 
This    is    Love's   equation    (see    [i;]).      Similarly,    we    can  derive   the 
boundary   condition 

^^°)  Vt   -  S*:cx  =    -P;'   Pxt 

which  is  to  be  satisfied  at  y  =  0.   Hero  p  =  p(0).   The  boundary 
condition  at  y  =  -h  becomes 

(11)  \i/  =  const  =  0  . 

The  boundary  condition  (10)  is  correct  only  if  P  and  vj;  have  the 
required  derivatives,  which  is  not  the  usual  situation  of  interest 
in  hydrodynamics.   It  will  be  replaced  when  necessary  later. 

4 

B.   VJe  will  now  investigate  the  possible  existence  of  waves  ^^Jith 
velocity  U,  i.e.  of  solutions,  \|/(x,y,t),  which  have  the  form 

^(x,y,t)  =  ei^^(^iUt)Y(y) 

With  k  real.   Such  waves  exist  if  one  can  find  non-trivial  Y(y) 
satisfying  the  boundary  value  problem 

(12)  (pY')'    -   k^pY   -   '^p'Y  =0,         -h  <  y   <  0 

(13)  Y'    -   [xY  =   0  at  Y  =   0 

(11;)  Y  =   0  at  Y  =   -h 

_2 

where   |x  =   gU      .      This   problem  can  be   regarded   in  two   ways:      either 

2 

\i   IS  held  fixed  and  X  =-lc   is  the  eigenvalue  parameter,  or  k  is 

held  fixed  and  |j,  is  the  eigenvalue  parameter. 


ry  n 


8    - 


WO.: 'Jo 
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First,  suppose  p,  is  held  fixed.  Vie   look  for  possible  values 
of  k  for  which  non-trivial  solutions  to  (12)  -  (ll|.)  exist,  i.e. 
for  possible  waves  having  a  fixed  velocity.   This  is  a  regular 

Sturm-Liouville  problem,  and  it  is  well  known  that  there  exists  a 

2 
discrete  spectrum  of  real  (simple)  eigenvalues,  ~^''-^^    '"^^^   associated 

eigenfunctions,  si  (y),  which  are  orthogonal  in  the  sense  that 

0 


y  p(y)f!'n(y)4^(y)dy  =  O   if  n  ^  m. 


-h 
The  eigenvalues  and  eigenfunctions  can  also  be  characterised  by 

means  of  a  variational  problem..   Let 

0 
(15)        (Y,L^Y)  =  r  (pY'^  +  l^p'Y^)dy  -  '^Q^Y^io) 

-h 
and  Q 


(16)       (Y,Y)  =  I!y|!^  =  y  pY^dy 


-h 

If  Y(y)  satisfies  the  boundary  condition  Y(-h)  =  0,  integration  by- 
parts  yields  _ 

(17)  (Y,L^Y)    =     '^      p(Y'2    -   2iiYY')dy 

^-h 
0 


^      p[(Y'    -   [iY)^   -   |iV]dy 


-h 
and   therefore, 

(18)  (Y,L^Y)>   -ix^IIyII  I 

The  eigenvalues,  -k  ,  and  the  eigenfunctions,  ^^,    can  be  obtained 
by  minimizing  (Y,L-j_Y)  with  \k\\^   =  1,  Y(-h)  =  0,  and  with  Y 
satisfying  the  orthogonality  conditions 

(Y,^^)  =  ...  =  (Y,ct^_-,)  =  0  . 


itcriJ- 
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Consequently,  there  exists  a  sequence  of  k   with 

IJ.  >  kj  >  k^  >  k2  >  ...  >  k^  >  ... 

p.  2 

and  ic  — ^  00  as  n  — >  co  .   Only  a  finite  number  of  the  k^  can  be 

n  ^^ 

positive,  and  therefore,  there  can  only  be  a  finite  number  of 
bounded  waves  (if  any)  for  each  velocity.   The  eigenfunction  4 

which  corresponds  to  the  shortest  wave,  has  no  zeros  in  -h  <  y  <  q; 
it  is  similar  to  the  usual  surface  wave  in  a  homogeneous  fluid. 
The  eigenfunction  ^     has  n  zeros,  and  at  least  n  internal  maxima 
and  minima,  and  the  corresponding  wave  can  be  called  an  internal 
wave. 

It  is  useful  to  examine  the  maximum -minimum  property  of  the 

2 

eigenvalues.   Every  eigenvalue,  -k^,,  can  be  found  by  taking  the 

least  upper  bound  of  the  minima  of  (Y,L,Y)  with  ||  y||  -j^  =  1, 

Y(-h)  =  0,  and  n  other  independent  linear  constraints  ([l]),ch.6), 

2 
Then  it  is  clear  from  (l5)  that  if  p'  <  0,  k^^  is  a  non-decreasing 

function  of  \x»      This  suggests  the  possibility  that  for  some  value 

2 

of  [J,  one  of  the  k   can  vanish;  such  values  of  \^   would  correspond 

to  critical  speeds.   This  will  be  clearer  when  we  examine  the 
problem  with  fixed  k  and  varying  p.. 

Remark.   Using  the  orthogonality  end  completeness  of  the  4^(y) 
makes  it  easy  to  write  down  a  proof,  similar  to  Weinstein's  for  a 
homogeneous  medium  [  9],  that  e"^  n^-^='^^"^Mj^(y)  (k^  real)  are  the 

ro 

only  bounded  waves  possible  among  -i}   with    \  p(y  )\l'(x,y,  t  )dy  <  qd  . 
One  need  only  expand  \1(  in  a  series  of  the  q  , 


.-;*      :  A-. 


*    .        » 
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*(?,y)  =  Tcn^^)^n^j) 


where  ^  =  x  ±  Ut,  and  the  c^  will  be  found  to  satisfy 

"    2 

c„  +  k  c  =  0  . 
n    n  n 

C.   Let  k  be  held  fixed,  and  let  [i   be  the  eigenvalue  parameter. 
The  problem 

LpY  =  -(pY')'+  k^pY  =  -ixp'Y 


with 


Y''  -  tiY  at  y  =  0   and  Y(-h)  =  0 


is  not  a  standard  Sturm-Liouville  problem  since  [i   occurs  in  the 
boundary  conditions.   However,  this  complication  introduces  no 
difficulties.   Let 

0 


(19) 

and 
(20) 


(Y,L2Y)   =     r      p(Y'^  +  kV)dj 


(Y,Y)   =    ||y|||  =    '      [  p'Y^dy  +    PqY^(O)    . 


(Here   we    require    p'    <  0. )      It   is    easy   to   see    that   for  Y  satisfying 
Y(-h)    =   0,    (YjL^Y)    is   positive  definite   with   respect   to    ||  YJp    . 
Integration  by  parts   gives 

IIyII  I  =   2     \     pYY'dy 


and   then 


_2  '^  0 

(21)        |(y12  ^    P(Y'^   +  A^)dy   >       ?    pY'2 

<-h  -Ih 


-2 


>  i      \'  pY'^dy 


^h 


dy 


0 
-^     ,.2 


V 


0       ->i 


2     \pYY'dy 

■<-h 


-h 


p 

^h 


1-1/2 


pY^dy       ?  pY'^dy 


(el) 


^u 


r  0    ~i,^/^ 


0 


1-1/2 


py  dy:  ■ 


-h 


J 


p(o)" 

.p(-iTi: 


1/2 


0 
Y'^dy 


-h 


1/2  '    0 


^   Y^dy 


J 


-h 


"•-1/2 

J 


rp(0)       il/2         7t 

^^r:ET^  2h 


One  can  conclude  (see  for  example  [  ^] )  that  there  exists  a  discrete 
spectrum  of  eigenvalues,  ^i,  =  lJ.j^(k  ),  with 


0  <  \i^  <  [i-^   <  [12 


IJt 


n 


and  \i     — ^   OD  as  n  — ^  oo  .   The  corresponding  sequence  of 
elgenfunctions,  Y  =  Y  (y;k),  are  complete  and  orthogonal  vjith 
respect  to  the  inner  product  (20),  i.e. 


<22).  (Y^,Y^)  =  - 


p'Y^Y^dy  +  PoY^(0)Y^(0)  =  0   if  n  /  m. 


-h 


An  arbitrary  function,    f(y),   defined   in   -h  <  y  <  0,    with    ||f  II2  <  00, 
can  be   expanded  in   a  series   with  respect   to   the   eigenfxmctions 
( for  any   re  al  k ) : 


CO 


(23). 

f(y) 

=  S       a      Y    (y) 
0 

where 

(21^). 

^n  = 

Y       -2 
^    2 

\    p'Y^fdy  +   p 

o  n 


All  these  conclusions  hold  for  any  fixed  real  k,  and  in  particular 
for  k  =  0.   This  proves  that  there  is  an  infinite  number  of 
critical  speeds,  U  =  [g/U„  (0)  j-'-^^,  and  U„  — >  0  as  n  — >  co  . 


n 


n 


n 


It  is  interesting  to  observe  some  of  the  features  of  (j.  as  a 

2 
f\anction  of  k   (which  characterises  the  dispersiveness  of  the 


.^  ■.., . 


s  te 


c-.-H 


,  3T  f/J.^.'   GwrTes   e  ni  L.  .■ 


i5    "T  --3   fv)*i 


t>f  Xfia'j   J30X.n  v,r;0  '-lol  tlod  a.no£Byionou 


a»j      ■;.* 


18     Ij 


10 


fluid)  without  making  any  further  restrictions  on  the  density 
distribution.   The  marimum-miniinuin  property  of  ^  ,  in  conjunction 

with  (19),  shows  that  M.  (k  )  is  monotonic  non-decreasing.   Ke  can 

2 
also  make  use  of  the  fact  that  ^x  is  an  analytic  fiuiction  of  1:  . 

A  formal  perturbation  then  yields 

0  0 

^'^^    dk^  =  1^  -  ?  -  ;?ii^ 

2  2 

which  shows  that  V-^i'k    )    is  strictly  monotonic  and  that  0   <  ]<.~   \i     <  c 

(c  >  0).   Furthermore,  since 

we  obtain 
(26)  3^  (-4)  <  0 


dk'-  k^ 

-2  2 

Consequently,  k   |jl  approaches  a  limiting  value  as  k  — >  oo  • 

2 
Formula  (25)  also  states  that  at  every  point  of  the    W-n^^"^  ^  curve 

the  slope  of  the  tangent  to  the  curve  is  less  than  the  slope  of 

the  line  from  the  origin  to  the  point. 

-1/2 
v^  -  Vj^(k;  -  gn^   ^  ^j^Q  progressive  waves  can  be  written 


in  the  form 


ik(x  ±   V  t) 

M'n  =  e        ^  y^(y;k) 


then  the   group   velocity/-   of  the  wave   i: 


(27)  Vg^=    (kv^).    =kv;.    v^ 


S-. 


11 


where  '  denotes  differentiation  '..-ith  respect  to  k;  (26)  implies 

that  V   >  0.   '^e  can  define  v  (k)  to  be  an  even  function  of  k, 
gn  n  ' 

monotonically  decreasing  from  U  =  v  (0)  to  zero  as  |k 1  — >  oo  . 

Consequently  kv  <  0,  and  it  follows  from  (27)  that  the  group 

velocity  is  always  less  than  the  phase  velocity. 

All  the  above  properties  of  u,  and  v  hold  for  all  n.   It 

n      n 

might  be  expected  that  in  some  iv-ays  the  surface  wave  will  differ 
from  the  internal  waves.   One  such  point  of  difference  is  the 
asymptotic  behaviour  (as  k  — >  oo  ) .   Since 

0  „0 

pS.   (Y'^  +  k2Y)dy  p(-h)  '^(Y'^  +  k2y)dy 


o 


(28)    _-J} <  —l^—  < 


(Y,L2Y)         Jj^ 


0  -  "iJYp"   ~       f^ 

r-L^Y^dy  +  PoY^(O)      ^        r2  \  Y^dy  +  P^y2(0) 

■Lh  -h 


where  r,  =  max  (-p')  and  r^  =  min  (-p')j  i^'©  can  obtain  bounds  for 

from  the  eigenvalues  of 
eigenvalues  of  the  problem 


LL  from  the  eigenvalues  of  two  auxiliary  problems.   If  a  are  the 
■^n  °  '^    '■  n 


(29)  Y"  -  k^Y  +  a  -1  Y  =  0 

(39')  Y(-h)  =0,   Y'  -  aY  =  0    Pt  y  =  0 


and  p        the   eigenvalues   of 


2„    .    .   ^2 


(30)  Y"    -  k^Y   +    PppHT^  "    ° 


Po 


(30')  Y(-h)   =   0,        Y'    -   P   pCrHy"^  =   °        at  Y  =   0 

then   °^n  ^  '-''n  —  ^n*   ■'"^  ^^  easy  to  obtain  the  asymptotic 
behaviour  of  a^  and  P^,  and  we  finally  get 


12 


2 
for  large  k  ,  and 

(31-)  0  <  A^  <  -^-^S-— ^  <  B„,    „  =  1,  2,  3,  ... 

2    .2  "^ 

for  large  k  +(^)  •   Thus  v^  is  of  the  order  of  the  square  root 

of  the  wavelength,  i,jhile  v  is  of  the  order  of  the  wavelength, 

as  the  wavelength  goes  to  zero. 


D«    It  was  shox'jn  in  the  previous  section  that  for  any  admissible 
p  (i.e.  p  >  0,  p'  <  0)  all  waves  travel  with  a  speed  less  than  U  , 
the  highest  critical  speed.   It  is  natural  to  ask  if  there  is  a 
wave  of  greatest  speed,  and  for  what  density  distribution  can  it 

occur.   This  is  equivalent  to  finding  [i" , 

[i'"  =  gib  [1^(0) 
among  all  admissible  p.   T-re  x,jill  show  that  |j,  =  h  ,  i.e.  that 
U  <  >/gh,  which  is  the  critical  speed  in  a  homogeneous  medium. 
This  can  be  easily  done  by  employing  Jacobi's  multiplicative 
variation. 

It  is  sufficient  to  prove  that 


(32)    Q(Y)  =  \   pY'^dy  +  h 


0  .0 

\   p'Y^dy  -  PqY^(O) 

-h 


>  0 


3X  ■  uy  -t-  r 
-h 

for  all  admissible  Y  and  p  (Y  is  admissible  if  It  has  a  piecewise 

continuous  derivative  and  if  Y(-h)  =0).   Let  the  function  Z(y) 

be  defined  by 

(33)  Z(y)  =  h(h  +  y)"^  Y(y) 

and  Z(-h)  =  hY' (-h) 

Then  Z(y)  is  continuous  and  has  a  piecewise  continuous  derivative. 

Substituting  (33)  in  (32)  gives 


rc'x 


'b£ 


to  > 


«».■■ 
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^k 


Q(Y)  =  N  p 
-h  •- 


(1  +  yA)^  Z'^+  ir^Z^  +  2h"^(l  +  y/h)ZZi 


0 


dy 


h"^  ]]  p'd  +  y/h)2z^dy  -  h-lp^z2(0) 


-h 

Integrating  by  parts  the  term  containing  ZZ',  x^e  obtain 
.0  0 


0 


2ZZ'p(l  +  yA)dy  =   p^Z^'CO)    -  \    p'(l  +  y/h)Z''dy   -  h"-^\    pZ'^dy 


-h 


-h 


-h 
and   therefore, 

0  ^  0 

(31;)  Q(Y)  =  W(l  +  y/h)^Z'2dy  +  h"^  ^  (y/h)  (1  +  y/h)p'Z^dy 

-h  -h 

Clearly,  Q.(Y)  >  0  since  p  >  0,  (1  +  y/h)  >  0  and  p'y  >  0,  and  Q, 

can  vanish  only  if  Z*  e  p'  =  0.   Thus,  Q  =  0  only  if  p  =  const 

and  Y  =  const  (1  +  y/h),  x-jhich  is  the  eigenf unction  for  the  problem 

(12)-(li|)  with  k  =  0  and  p  =  const.   Therefore, 

gib  ti^(O)  =  [i''"  =  h-1  or  lub  U^  =  U'"'  =  ,/gH  . 
P  P 

It  might  be  noted  that  as  h  — >cei  ,  \i"   — >  0  and  U'"'  — >  oo  ,  i.e. 
there  is  no  critical  speed  for  an  ocean  of  infinite  depth  and  p  = 
const. 


E.    For  an  ocean  of  infinite  depth  the  situation  is  much  more 

complicated  since  the  Sturm-Llouville  problem  (12)-(li;)  is  singular. 

V/e  vjill  not  discuss  this  problem  here,  except  to  give  an  example 

with  continuous  p  which  has  an  infinite  number  of  critical  speeds. 

This  may  have  some  importance  since  it  suggests  the  possibility 

of  the  existence  of  solitary  waves  in  an  infinite  domain. 

jp  e~Py,  -  h  <  y  <  0 
Let  p(y)  =1    _g 

'  p  e  ^    ,    -  oo<  y  <  -h 


d-': 


■■{^' 
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Below  y  =   -h,    Y  =  Y-,    satisfies 

(3U)  YiJ   -  ]^h^  =   0 

XTihile   for   -h  <y  <  0,    it    satisfies 

(35)  Y"    -   PY'    +   p-PY  -  A  =   0    . 

Therefore, 

Yi   =   ell^ly 

(in  fact  Y^   -    |k|Y  =   0     for     y  <  -h)    . 

In   -h  <  y  <0, 

Y  =   e-*-/^   P^[A   cosh   a(y  +  h)    +  p   sinh   a(y  +  h)] 


where 


2        1      P  2 

a     =  ^  P'^^   -  lip  +  k"^    9 


Since  Y  and  Y  must  be  continuous  at  y  =  -h,  we  can  evaluate 
A  and  B.   The  free  boundary  condition  then  yields  the  following 
equation  to  determine  \i: 


(36) 


tanh   ah  _  ______  ji.  -    |k  | 


"aET 


|kh|(  Ik  I   +  ^  p-tj.)    -  ^  tiPh 


The   critical   speeds   are   obtained  from  this   equation  with  k  =   0; 

(37)  tanh  ah  =   -  2p"-^a 

2        12  . 

and  a     =   j-  p      -   [iP    . 

There   are   an  infinite  number  of  solutions,    [o.   (0),    with  all   the 
1X^(0)    <l  p. 
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Section  II     Travelling  surface  pressure;  steady  solution. 

A,    In  this  section  we  will  consider  the  disturbance  created  by 
a  moving  pressure  distribution  applied  at  the  surface.   Suppose  a 
steady  pressure  is  applied  at  t  =  0,  and  moves  with  a  constant 
velocity,  U,  in  the  direction  of  the  negative  x-axis,  i.e. 

(38)  P  =  P(x  +  Ut)  =  P(^),   5  =  X  +  Ut. 

We  assume  P  is  continuous,  and  vanishes  outside  of  the  finite 
interval,  |^|  <Jl,      The  derivative,  P'(^),  need  not  be  continuous; 
however,  we  will  require  P(^)  to  behave  no  worse  than  ij?-^^)^, 
0  <  e  <  1,  as  f;  approaches  either  end  point  from  the  interior. 

It  is  natural  to  expect  that  when  viewed  from  a  coordinate 
system  moving  with  the  pressure,  the  disturbance  will  tend  to  a 
steady  state.   Consequently,  we  look  for  a  solution  of  the  form 
'if   =^(?,y).   Introducing  this  into  (9)  and  (10)  of  the  preceding 
section  yields 

(39)  u2tpA¥+       P'lylfJ        .      gP%£       ,        0 

The  BC  (1|0)  is  not  suitable,  since  even  P^.  may  become  infinite 
at  ^  =  +i/.   However,  we  can  write  it  in  the  form 

%  -  ^^*   P^  'h,   =  0 

and  try  to  satisfy  the  condition  that  the  quantity  inside  the 
brackets  vanish  (If  we  go  back  to  the  original  equations  and  look 


r-    i.-aoii-:;U'd.      *'- 


,-kn    ; 


<'•  -f  1^  *j.X  i      l^i- 


:h)-0 


~  i-C*" 


1      •• 


t'lol-^'i^^ai 


'^rt     r-.-.-  -!     >■ 


■^.c*ii:.^c>etc    t::i;t    V 


•ad-lun^. 


«   \h 


1? 


for  a  solution  which  depends  on  f"  and  y  only,  the  superfluous 
differentiations  will  also  not  appear).   We  obtain  the  problem 

(III)  (pfC-)^   +   (Pty)y   -   IXP'T=   0,  tX   =   g/U^ 

ih2)  Uy  -  [i^^=   -  p^  at  y  =  0 

(U)  "'i"=  0  at  y  =  -h. 

As  is  common  in  situations  of  this  type,  the  problem  defined 
by  ( i^.1 ) - ( Ij.3 )  has  in  general  more  than  one  solution.   It  is 
clear  that  for  speeds  below  the  highest  critical  speed  (for  the 
given  p),  all  the  possible  waves  are  solutions  to  the  homogeneous 
problem  (P  =  0),  and  can  be  added  to"!".   To  insure  uniqueness 
one  can  introduce  an  artificial  "viscosity"  and  find  the  limit 
solution  as  this  parameter  vanishes  (Rayleigh),  or  one  can 
impose  conditions  on  the  behaviour  of  the  disturbance  at  oo  . 
We  will  require  thatV"vanish  far  ahead  of  the  disturbance 
(^  —>  -00  ) ,  and  that  it  be  bounded  as  E,   — > +oo  .   Then,  uniqueness 

r^        .-2 
among  the  class  of  functions  with  \  p{j)^[j   (€>7)dy  <oo  ,  for  all  ^, 

-h 
can  be  easily  proved  by  using  1-reinstein's  procedure  (see  remark 

in  IB).   A  more  satisfactory  justification  can  be  made  by 

considering  the  time  dependent  solution,  and  showing  that  it 

approaches  the  steady  solution  (see  III). 

B.    The  solution  to  the  problem  (ii.l)-(i4-3)  with  the  condition 
that'll  vanish  far  "upstream"  can  be  obtained  by  a  superposition 
of  solutions  to  (Li-1)  which  are  harmonic  in  £  and  satisfy  the 
boundary  condition  at  the  rigid  bottom.   We  know  from  previous 
results  that  e"""  -Y(y;k)  is  a  solution  to  (l|l)  if  it  satisfies 


..  ^■-■P ; 
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ikk)  (pY')'  -  k^pY  -  i^p'Y  =  0. 

Choose  Y  so  that 

ik^)  Y  =  0   and  Y'  =  1    at  y  =  -h. 

Then  Y  is  uniquely  determined.   It  should  be  noted  for  later  use 

that  Y(y;k)  is  an  analytic  entire  function  of  k  for  each  y.   Now 

consider  a  superposition  of  such  solutions  for  all  possible  k: 

T(^»y)  =  5A(k)e^'^^%(y;k)dk 
L 
where  the  path  of  integration  in  the  complex  k-plane  will  be 

fixed  later.   To  satisfy  the  boundary  condition  at  y  =  0  formally 

■\]r^(?,0)  -  |aT(^,0)  =  ^A(k)[Y'  -  tiY]^  e^^^dk  =  -  -^  ?{E) 

Jb 

where  the  subscript  o  doubles  the  value  at  y  =  0,   If  the  path 

is  such  that  the  Fourier  Integral  Theorem  can  be  applied  (formally) 

we  can  evaluate  A(k),  obtaining 

(i;6)    T(e,7)  =  -  (2^PoU)-^  Y[V^   Vj^^    e^^'^ik 

L  ° 

where 


P(k)  =  \p(r)e"^^^dE  . 


.1 
The  denominator  of  the  Integrand  has  an  infinite  number  of  zeros 

on  the  imaginary  axis,  and  if  [i  >   ^  (0),  a  finite  number  on  the 

real  axis.   Suppose  ^j,  >  1x^(0)  but  is  not  one  of  the  critical 

speeds,  and  let  k  =  +  k  ,  n  =  0,  1,  2,  . . . ,  N  be  the  zeros  on 

the  real  axis  (k^^  j^   0).   Then  for  L  we  pick  the  path  which 

coincides  with  the  real   k-axis  except  near  k  =  +  k  ,  where  it 

consists  of  semicircles  in  the  lower  half  plane.   It  will  be 
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proved  now  that  T."",  as  given 
by  (L|-6)  is  indeed  a  solution 
to  ( 1^.1 )  -  ( [j.3  )  and  satisfies 
the  required  condition  as 
4  — >  -CD  . 


-^o  -^1 
L 


(k) 


^1   ^o 


C.   To  prove  that  ([|-6)  is  a  solution  we  need  the  behaviour  of 

the  integrand  for  large  |k:|,  k  real.   The  asymptotic  behaviour 

-1/2 

of  Y  is   easy  to   obtain  if  we    change  variables.      Let  V  =   p     '    Yt 

Then    ik-k)    s^d    (L|.5)   becoiiie 


ikl) 


V"    -    l^  £V  ^  ^'-  lif]^  H.   k2]V  =  0 


ik^)  V(-h)  =  0,  Vi(-h)  =yp(-h)  =  1  if  we  normalize  p. 

Then  it  is  easy  to  see  that 

(L^9)        Y  =  p"^Ar^  sinh  k(y  +  h)[l  +  0(k~^)] 

Y'  =  p"-'-/^^  cosh  k(y  +  h)[l  +  0(k~^)] 
and 

Y'  -  nY  =  p"-"-/^  cosh  k(y  +  h)[l  +  0(k"^)] 

uniformly  for  -h  <  y  <  0  and  large  |k|.  Therefore, 


(50) 


Yfvk')        ^O  pl^^'y  -l^,' 


From  fie  assumptions  made  about  P(^),  we  know  that  P(k)  = 
0(k"    ),  0  <  e  <  1,   Consequently,  the  integrand  is 
exponentially  bounded  for  large  |k|  and  y  <  0,  and  therefore, 
the  integral  in  (1|.6)  is  absolutely  convergent.   Furthermore,  it 
can  be  differentiated  any  number  of  times  x^>ith  respect  to  t, 
and  y  under  the  integral  sign,  and  this  proves  that  (I1.6)  is  a 
solution  to  the  differential  equation  (l+l). 
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To   shoi-j  that    (1;6)    satisfies    the   boiondarAr   condition   (i|2), 
form^     -   ni^-^for  y  <  0,    and  let   y  — >    0. 

V 


^^-^v=-^S^w^''''^^^%- 


Li 


.-l-e^  ..^  Y'  -  H-Y 


Since  P(k)  =  0(1^""  "  )  and  ryT—Yl^  certainly  stays  bounded  as 

^  -'o 

y  — i>  0,  the  integral  on  the  right  hand  side  is  absolutely 

convergent  for  y  =  0,  and  therefore  we  can  go  to  the  limit  under 

the  integral,  obtaining 


lira   [fy   -   ^f]    =   ^  ^  jF(k)eil<^  dk. 


y=0      '-^  -'^^^o      L 

F(k)    is    analytic,    of   course,    and  we    can  deform  L  to   the    real   axis 


CO 


P(£) 


e^l^?  dk 


-00 

1 

PoU 


which  proves  that  (i;2)  is  satisfied.   Obviously  (I|-3)  is  also 

satisfied. 

Before  going  on  to  the  asymptotic  behaviour  ofljj{C,fj),    we 

note  that  for  y  <  0  it  is  certainly  permissible  to  substitute  the 

expression  for  P(k)  into  ([|.6)  and  then  to  change  the  order  of 

integration,  which  results  in 

i 


(52)   "5r(£,y)  =   -  ~jj   ^P(s)  K{£,s;y)ds 

-i 


(53)     K(?,s;y)  =  ^i^ll^jei^^^-^^dk 

is  the  stream  function  due  to  a  6-function  in  the  surface  pressure. 
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In  the  neighborhood  of  s  =  4  ^-^d  y  =  0,  k  has  a  logarithmic 
singularity. 

It  will  now  be  shown  that  'vj/vanishes  as  E,   — >  -co.   Let  the 
integral  in  ([|.6)  be  expressed  as  the  sum  of  two  Integrals,  one 
over  a  portion  of  the  real  k-axls  and  the  other  over  the  2(N+l)  semi- 
circles.  Since  the  integrand  is  0(k~  ~  )  for  large  |k|  on  the 
real  axis,  the  first  contribution  vanishes  as  |^|  — >oo  by  the 
Rieraann-Lebesgue  Lemma,   To  evaluate  the  second  contribution, 

consider  one  of  the  semicircles,  say  the  one  near  k  =  k  ,   The 

'  n 

integrand  can  be  written  in  the  form  f (yjk)e        ,  with  f  bounded 

on  the  semicircle.   On  the  semicircle,  k-k  =  a(cos  ©+i  sin  e), 

'    n  ' 

a  >  0,  TC  <  6  <  271,  and  so  we  obtain  an  integral  of  the  form 

2k 

fi(y,e)e-°^  ^^^  ®-«  d9 

with  f,  bounded.   As  4  -s-  -oo  the  integral  vanishes  like  IcI"   . 
Since  there  is  a  finite  number  of  semicircles,  the  second  contribu- 
tion also  vanishes, 

D,   To  complete  the  proof  that  (l|-6)  is  the  steady  state  stream 
f\inction,  we  now  want  to  examine  the  principal  effect  of  the 
surface  disturbance:   the  waves  created  in  the  wake  of  the 
disturbance.   The  surface  elevation,  \"i    is  given  by 

yi  =  -  i  lim  '|-(>?,y)  . 
^       y=0 


I 


•J--V- 
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Consequently  we  will  examine  ^<  as  <^  — >  +00  ,   In  order  to  do  this 
we  need  the  behaviour  of  [Y'  -  |iY]   in  the  complex  k-plane  near 

its  zeros,  k  =  i  k  ,  and  this  we  proceed  to  do  first, 

2 

Suppose  9{j)   is  an  eigenfunc tion  and  -k  <  0  is  the  corres- 
ponding eigenvalue  of 

i^k)  4^  =  -(pY')'  +  Hp'Y  =  -k^pY 

(55)  Y(-h)  =  [y'-p.Y]^  =  0  , 
i.e. 

(56)  44.  =  -k^ocj, 

(57)  4(-h)  =  [ct)'-^(l)]o  =  0 

and  normalize  (j)  to  have  (j)'(-h)  =  1,   Now  let  k  be  a  complex 

2 
number  in  the  neighborhood  of  k.,,.,  and  let  Y(y;k)  be  the  solution 

of  (54)  which  satisfies 

Y(-h;k)  =  0  ,   Y'(-h;k)  =  1  . 
The  differential  equation  {Sk)    can  be  written  in  the  form 

(58)  L^Y  +k|pY  =  -(k^  -l4^>^  • 
Consequently,  (Sk)    can  be  replaced  by  the  Integral  equation 


(59)   Y  =  c|,(y)  -(k^-k|)  r  [t(y)X(»l)  -<l)(v;)X(y)lp('))Y('{,-k)dv^ 

-h 

where  ^(y)  is  the  solution  to 


f  ' » 


■f  ,v !  I 


{  'r 


,M-" 
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2 


(61)  Xi'h)   =  -  -^X^   ,   Z'(-h)  =  0  . 


Then 


(60)  L^X   =  -  K^X 

y  ^ 

(62)  y'  -ixY  =  (cj)'  -n(j))  -(k^-k|)  r  [[<!)' (y)  -txc|>(y)]X(»^) 

-h 

-  [^'(y)  -txX(y)](f(Vi)yi*(r|)Y(y|;k)dVi  . 
Letting  y  — >  0,  we  obtain 

(63)  lini  [Y'-tiY]  =  [Y'-jiY]^=  [(j)-tx(})]^ 
y=0  0         0 

0 

-(k^-k?)  \    l^'  -ii^]^XiV()-[K^  -[i^]JiY^)^{^n{vi;^)dYi 
-h 

0 
=  (k2-k|)[;^'  -(i^]^  V  0(»])cl)(v|)Y(i^|k)dy| 

[X'  -[iX]     cannot  vanish,  and  Y(y;k)  — >  <})(y)  uniformly  in  y  as 

2  2 

k     — >kTI.      Therefore, 

0 

(6if)      [Y'-txY]^=    (k2-k^.)[  X'-p,;^]^    j    p(>j)(l)^(v^dv^+   0((k-k.^)^) 

-h 

since  Y  and  Y  are  analytic  In  k  for  every  fixed  y.   This  shows 
that  [y'  -|j,Y]  has  only  simple  zeros  on  the  real  k-axis  (provided 
\i   does  not  correspond  to  a  critical  speed).   We  can  also  evaluate 
[X  -[iX]      since  we  know  the  values  of  the  V/ronskian,  W: 

W(y)  =  tX'  -//  =^  . 
Therefore, 


)     f  .■         •       \j(.^ 


y  /    ■■■■ 


rr- 


■'G    Vj 


*  •  *     *  ■ 


-■  \  "-; 


2k 


W(0)  =  [X'  -(xXJot^o)  =  J- 


and  this  yields  finally 


(65) 


(k^-k!)   0 


[Y'-tiY]   =— r-rro?-  j  P{y)<t'^  (y)dy  +  0(  (k  -  k  J^)   . 

'"°^     -h 


E,   In  order  to  evaluate  lira  "'i  ('?,y),  we  can  deform  the  path  L 

C->+oo 
+ 

into  L  ,  which  coincides  with  L  on  the  real  axis,  but  passes  above 
the  singularities.  Clearly, 


L   L 


+  271 


residues  at  k  =  k. 


n  • 


As  4  — ■>+oo  the  integrand  is 
exponentially  bounded  on  the 
semicircle  of  L  ^  and  con- 
sequently I  — >  0,  leaving 

only  the  contribution  at 
the  singularities.  Using 
(65)  we  can  compute  the 
residues,  obtaining  (as 

€  — >  +CO  ) , 


(k) 


T       Y(Ojk^)Y(y;k^)  IkC 

(66)    l/(C,y)  ^>  ^£Z "0^^ P(kn^^ 


0 

f 

-h 


^  K  I  >'Y^(7;k^)dy 


)  + 


,  •    '  ■-. .  t.' 
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where  the  sxaranatlon  extended  over  all  the  poles  (k  =  i  k  . 

^  =  0,1,,.,, N),   Since  Y  Is  an  even  function  of  k,  we  can  combine 

the  contributions  at  +k^  and  -k^: 

p 
T      Y(0,-k  )Y(yjk  )     ^ 
(67)   >i'(€,y)  -»U"^  ^  ^ 2 j  p(s)  sin  k^(C-s)ds 

^n  {  P^^^1'>\,)^1   "^ 

where  the  svmimation  is  now  over  the  positive  k  ,  or  using  the 
designation  ^^il)   for  Y(7;k^)  as  in  IB, 


<t.^(0)<j,^(y) 


^ 


ds 


(68)  f(C,y)  -^U-^EZ  ^  „  ;"  .   f  P(s)  sin  k^(?-s) 

^n"^n"   ^JL 
and 

?      ^ 

(69)  ^(€)  -^  -  -^EZ  --TT-rr^  i  P<^)  si^  k  (4-s)ds 

as  C  — 5»  +00  • 

We  now  have  a  picture  of  the  distiorbance.   When  the  pressure, 

P,  moves  at  a  speed  greater  than  the  highest  critical  speed,  U^ 

(l,e,  n  <  M-  (0)),  there  are  no  waves  far  ahead  or  far  behind  P, 

If  U,  <  U  <  U  ,  there  is  one  wave:  its  length  decreases  as  U 
1    ,    o 

decreases.   Below  U,  a  second  wave  appears;  its  length  is  greater 
than  that  of  the  first  wave,  and  also  decreases  as  U  decreases. 
Below  Up  a  third  wave  appears,  and  so  on.   Par  away  these  waves 
behave  like  the  progressive  waves  corresponding  to  the  speed  U, 
It  should  be  noted  that  the  wavelengths  are  not  harmonically 


•  » 


•Xv 


v:>*X.I 


'(^      P,  - 


26 

related  in  general,  and  the  disturbance  behind  P  is  therefore  not 
periodic  in  ^» 

P.   The  waves  behind  the  surface  pressure  carry  off  energy, 
giving  rise  to  a  "resistance",  R,  which  must  be  overcome  in  moving 
the  pressure  distribution.   In  the  case  of  constant  density  the 
resistance  is  proportional  to  the  height  of  the  (single)  wave 
train  far  behind  the  disturbance.   In  the  present  case  the  internal 
waves  also  contribute  to  R,  and  their  effect  can  be  sizable  as 
noted  in  the  Introduction, 

The  resistance  can  be  foimd  by  summing  the  horizontal 
components  of  the  pressure,  P,  over  the  surface  deformation,)^, 
created  by  P,   In  the  linearized  approximation  this  gives 

(70)  R  =  -  <:  P(?)>('(?)d? 

and  in  view  of 


.f 


.-1^^    _x  _  -,,.„    1 


(71)  yAa   =  -  lim  U"'-ar(^,y)  =  lim  —^    \    P(  s  )K(4,s;y  )ds 

^        y=0  y=0  p^u2  J 

we  get  J  -a 


(72)  R  =  - 

ji 


^  lim  rP(?)dS  V  P(s)K^(?,s;y)ds 


Ve  cannot  set  y  =  0  linder  the  integral  sign  because  of  the 
singularity  in  K  (^,s;y),  which  behaves  like  a  double  layer  charge 
distribution.   However,  the  limit  can  be  easily  obtained.   Let  K_ 
be  split  up  into  an  even  and  an  odd  function: 
(73)  K^(^,s;y)  =  E(£,s;y)  +  0(^,s;y) 

where         E(^,s;y)  =  ^[K^(^,s;y)  +  K^(-?,-s;y)] 

and  0(?,s;y)  =  |[K^(^,s;y)  -  K^(-?,-s;y)] 


if.'-'-' '    ■• 


.>-■ 


lio  wc 


U  -   -   (^)v( 


« 


^t<'K^  I  -J-" 


'CS) 


iv;^-,^.-)^}t  ~  (y;j^.;^)3^U|  =  {y.iB,,'i)o 


Zl 


Clearly,    for  every  y  <  0  the   odd   pi  rt   of  K     does  not   contribute 

to   the    integral;    thus  n      j\ 

(7i|-)  R= i^  lim  S     ^     P(?)    P(s)   E(^,s;y)d.^ds 

making  use   of   (53),    we   obtain 

(75)  E(4,s;y)  =  ^^  ^^^7^  °°^  ^'^^"•^^  ^^ 

L  ° 

Since  Y  and  Y'  are  even  functions  of  K,  the  integrand  is  odd, 
and  the  integral  over  the  part  of  L  vjhich  coincides  with  the  real 
axis  vanishes.  What  remains  is  the  integral  over  the  semicircles, 
which  can  be  evaluated  since  we  Icnow  the  residues  at  k  =  +  k  . 
Carrying  out  this  computation  results  in 

-r —     <l>    (0)    (|)    (y) 

(76)  E(?,s;y)   =  -  1  p^    >_       ^ '^       g cos  k^(?-s) 


and 
(77) 


_  c  c  T —       ^  (0)4*  ( y ) 

R  =  ^  llm     \  \     P(?)P(s)  2_  "    „,   "g  -  COS  k   (?-3) 


-^Z  \S     P(.)    P(s)[Reei^(^-)^,i^ 

4!(o) 


2U^  Fr'>   0    II 


|P(k    )l 


2 


>    I )     1 1   ai      ii_ 


where,  as  before,  P(k)  is  the  Fourier  Transform  of  P(^)  and  the 


summation  is  over  the  positive  k  . 
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G.    It  is  natural  to  ask  whether  R  can  vanish.   For  the  case  of 
constant  density  one  can  argue  this  possibility  verj'-  siraply.   Take 
any  tijo  identical  pressure  distributions  one  half  of  a  wavelength 
apart.   The  trailing  waves  produced  by  each  one  will  differ  in 
phase  by  l8o   and  will  cancel  each  other  far  downstream. 
Consequently,  the  waves  '.'ill  vanish  both  upstream  and  dovjnstream, 
and  the  resistance  must  vanish,   'Tien  there  are  two  or  more  waves, 
and  the  resulting  wave  train  is  not  periodic,  this  argument  fails. 
However,  (77)  shows  clearly  that  R  =  0,  if  and  only  if,  P(k  )  =  0 
for  n  =  0,  1,  ..,,  N.   This  In  effect  constitutes  2(N+1) 
orthogonality  conditions, 

P(5)[sin]  k„^.  d?  =  0,     n  =  0,  1,  .,.,  N. 


5, 


n 

There  is  an  infinite  number  of  independent  dif ferentiable  functions 
in  the  space  orthogonal  to  cos  k^  and  sin  k  ^  (n  =  0,  1,  ..,,  N) 
over  (-",'<-),   A  linear  combination  of  any  three  of  them  can  be  made 
to  vanish  at  f!  =  +  X,   Consequently  there  exist  pressure 
distributions  of  the  type  we  have  admitted  which  produce  no 
resistance  (at  one  speed). 


Remark:  The  resistance  will  never  vanish  (below  the  highest 
critical  speed)  if  P(k)  >  0  for  all  k.   One  can  construct  examples 
for  which  this  is  true.   Take  any  continuous  and  dif ferentiable 
Pi(^),  which  vanishes  outside  of  |f  |  <  iJ-,   Form  the  convalution 


\   '"'  '^^  ~        V^(^)  ^^^  -  s)  ds    and  then  square  it.   The  function 

ca] 


(Pjl^-'J-  Pj^)   is  an  admissible  pressure,  and,  as  can  be  easily  computed. 


its  Fourier  Transform  is  positive  for  all  k, 


(•»5-)   This  example  was  suggested  by  Dr.  R.  Silverman, 


>XdIi 


.a  , 


29 

Section  III,   Travelling  surface  pressure;  time  dependent  solution, 

A,   Suppose  the  fluid  is  at  rest  and  that  a  travelling  surface 
pressvire,  P  =  P(x+Ut),  is  applied  at  t  =  0,   We  want  to  determine 
the  time  dependent  stream  function,  i.e.  determine  ^(x,y,t) 
satisfying  the  differential  equation  (9)  and  the  boundary  condition 
(10),  (11),  the  condition  that  ^  be  boxmded  at  oo  ,  and  the  Initial 
conditions 

(78)  ^(x,y,0)  =  \l'^(x,y,0)  =  0  . 

When  P  is  not  sufficiently  differentiable,  (1)  is  replaced  by 
another  condition,  which  need  not  be  specified  for  the  present 
problem).   It  is  clear  on  physical  grounds  that  both  of  the  condi- 
tions (78)  are  required,  since  the  initial  density  distribution 
must  be  known  as  well  as  the  initial  velocity  field. 

Since  the  steady  state  stream  fimction,  ^^(x+Ut,y),  given  by 
(1|.6)  or  (50),  satisfies  the  differential  equation  and  all  the 
boundary  conditions,  it  is  necessary  only  to  find  a  solution  \I/'  to 
the  problem  with  the  homogeneous  boiindary  condition 

(79)  il/;^^  -  S^ic  =  °   at  y  =  0 
and  the  initial  conditions 

(80)  ^  \J/'(x,y,0)  +  ^(x,y)  =  0 

\I;^(x,y,0)  +  Ui?'^(x,y)  =  0  , 
Then, 

ilf(x,y,t)  =  \l/^x,y,t)  +\r(x+Ut,y)  . 

The  fixnction  v]/     can  be  found  by  superposing  known  solutions  with 

2  ikx  -^^^   * 

V     =  g/|J.   ,   we  know  that   e        e  ^  Y   (y;k)   is   a  solution  of  the 
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differential  equation  with  the  homogeneous  boundary  conditions 
provided  Y  (y;k)  satisfies 

(81)  (pY^)  ♦  -  k^^Y^  =  n^  p'y^  =  -^  p  Y 


^n 

(82)  Y^  -  ^  Y^  =  0         at  y  =  0 

^n 

(83)  \=  ^  at  y  =  -h  . 

Now,  It  can  be  easily  verified  that  a  solution  to  the  problem  is 
given  (at  least  formally)  by 

(8^)     ^  =  -(2.^^U)-1  J  m^^f^  e^'^^-^^^^    dk 

L         o 
ikx 


with 


T   r        ii^v^t  -ikv  t  ^  ■ 


n 

5^(k,0) 
Since  ^  (k,0)  =  1  and  — S-_ =  ikU,  the  initial  conditions  will  be 

satisfied  provided 

00 

(85)  Y(y;k)  =  >___  a  (k)Y^(y3k) 

0 

where  the  Fourier  coefficients,  a  (k),  are  obtained  from  (23), 

In  view  of  the  definition  of  Y(y;k)  as  the  solution  to  {ki\.) 

taking  on  the  initial  values  ik-S) »   we  can  evaluate  the  a  (k)  by 

using  Green's  Formula: 

0  0 


^  (Y.V-YL^Yjdy  =  [pCYY^-Y'Yj]  j 


0   -h 


Therefore  "" 


J  P'yy^dyf  ^^^Y(0)Y^(0)  =ej^(0)         <^ 
-h  ^ 


I    'i    yv- 


0    ^   -i    .SS.- 


el  m^L^xstq  crl^'t  o'   ,<•.  5-%-.'^, 
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and 

T  p  Y  (0)[Yt-^xy]  tY^(0)[Y'-p.Y] 

(86)  a^(k)  =  (,^-  ^)-l  ^ ii-^  =  (,^.,)-l_^  o 

-yY^dy  ^  p^Y^(O) 
-h 

Thus  we  finally  obtain 

(87)  \]/(x,y;t)  =F(x  +  Ut,y)  +  ^^^{x,j;t) 


dk 


(87-  )    ^(x  .  Ut.y)  =  -  (2,p„U)-^  5  Z||llto|^^!!±:!!!' 

if     ik^  V^  Y„(0;k)Y  (y;k) 
(87")    >lr'{x,y;t)  =  (2nU)-M  FCkje^*"^  2_  ^^ "   „g^(k,t)dk 

This  is,  of  course,  only  a  formal  solution. 

It  is  interesting,  at  this  point,  to  compare  (8?")  vjith  the 
analogous  result  for  the  constant  density  case.   When  p  =  const. 
\}('  satisfies  the  potential  equation  and  only  the  first  term  in  the 
sum  in  (67")  remains;  [x  is  determined  from  [i   tanh  kh  =  k  and  Y  = 
sinh  k(y  +  h).   The  resulting  integrand  contains  an  exponential 
factor,  e'  '•^,  which  permits  one  to  differentiate  the  integral  any 
number  of  times,  and  this  makes  it  easy  to  prove  that  it  represents 
a  "genuine"  solution  to  the  problem.   For  p  r   const  the  situation 
is  more  complicated.   It  is  difficult  to  estimate  the  asymptotic 
behaviour  of  Y   (i.e,  for  large  n  and  all  k),  and  the  Y^  are  not 
exponentially  damped  in  any  case  (n  >  1).   This  brings  up  the  usual 
difficulty  in  the  method  of  separation  of  variables  for  time 
dependent  problems.   One  can  prove  that  the  right  hand  side  of  (87") 
is  differentiable  a  sufficient  number  of  times  by  requiring  that 
P(k)  should  go  to  zero  sufficiently  fast  as  |k|  — >  co  ,  which 
imposes  smoothness  conditions  on  P(^)  and  its  derivatives  in  excess 
of  what  is  reasonable  for  the  hydrodynamical  problem. 


ve  will  consider  only  the  special  case  where  p  =  p^e 
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(p  =  const).   This  reduces  ( 8l )  to  a  differential  equation  with 
constant  coefficients  and  makes  it  easy  to  obtain  the  asymptotic 
behaviour  of  Y  .   In  the  next  section  (B)  it  will  be  shovm  that 
the  right  hand  side  of  (87")  represents  a  continuous  function, 
differentiable  with  respect  to  t  and  satisfying  the  initial 
conditions  (80),  and  in  section  C  that  \]/ '  -^  0  for  every  fixed  y 
and  X  +  Ut  as  t  — >  00  . 


B.     If  p  =  p  e"^^,  (81)  reduces  to 
(88)     Y"  -  pY'  -  k^Y  +  PtiY  =  0 

The  function  e"^   sinh  p(y  +  h),  with  a  =  (jj|{3  +  k  -  pix)  ^  , 
satisfies  (88)  and  the  condition  Y(-h;k)  =  0.   The  free  surface 

condition  then  yields 
/Ao\      tanh  ah  _      1 

which  determines  p^  as  a  function  of  k  .  l-'e  will  need  the 
behaviour  of  ^^  and  Y^^  when  \i^   is  large,  i.e.  for  large  k  or  for 
large  n.   Prom  (39)  it  is  easy  to  obtain  the  following  approximations 

(90)  ^^   =  Iki  +  0(e-2|^hl) 

1 

y^  =  e-2^^  sinh  |k|(y  +  h)  +  0(e-l^^'^ 

for  large  [k| 

(91)  \i^   =  P"^[k2  +  (^)^]  +  0(1) 

ipy  -, 

Y^  =  e'^    sinh  a(l  +  y/h)  +  O(i-) 

for  large  n,  uniformly  in  kj  n  =  1,  2,  , e .  » 

Then, 

Y  (0;k)  Y  (y;k)       k  1 


•    .    «    -1. 


=:         Y 
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Y^(0\]^)  Y^(y;k) 


—  Of    •»  \n^     n-K  sin  mi(l  •<-  y/h) 


Let   ^JL  be   fixed.      Then,    as  was   shown  in   section  1,    there   is 

some   N   such  that   |i   -   [i  >  0  for  n   >  N   and   all   real  k.      Consequently, 

(87")    can  be   rewritten  in  the   form 

/Mr  00  00.  "^) 

(93)  vlf'    =    (2t:U)'^|    Z_J(    )dk  +^      '^^  ^  ^    ^' 

^      L  -^-00 

Since  each  term  of  the  integrand  is  analytic  in  k  in  the  neighbor- 
hood of  the  real  axis.   Prom  (8I4.)  it  follows  that  for  real  k 

|p^(k;t)|  <  1  +  JI  =  1  +  U(tJ.^/g)^/2 
Consequently, 

Y„(0;k)  Y„(y;k)  p^^kjt)e^  /       . 

CO 

for  sufficiently  large  n.   Therefore,  2   A   can  be  made  arbitrarily 

M   ^ 
small  by  choosing  M  sufficiently  large.   Furthermore,  since 

GO  CO 

P(k)  =  0(k"""''"^),  \        F(k)(^^  Ara^^^  °^^  ^®  "^^^^  arbitrarily  small 

-co 

by  taking  M  sufficiently  large,  which  shows  that  the  integration 

and  summation  in  (87")  can  be  interchanged: 

^  ^-"1  r  00  00       f 

Ok)         r  =  (27iu)"^\  H  J  (  )<ik  +  n  \     (  )dk  4 

Since  \5;'  is  the  uniform  limit  of  continuous  functions,  it  is 
continuous  in  (x,y,t);  this  shoxx'S  that  the  first  of  the  initial 
conditions  in  (80)  is  satisfied. 

To  prove  that  >!/'  can  be  differentiated  with  respect  to  t,  we 
observe  that 


n 


■en.:;     ^e-io'-rnscii'Zif'i      .  ^t^o 


•.--  -,  r    ..  r 


<» 


?.< 


CO- 


.1    cv.vmo   doiriw    ^t:^ri'^:j::X  y^^^^^-'^J^^'^^J-'B 


00      '.:/j 


cx>- 


■(uttS)  «  »<r 


..   -ini/l  ^isotsn.t^fnofi  lo  tXixilf.  nr-  "  '     •  sr:. 


8m    ,J  o;^    ■  --^t  rf.tlw  food'B.tJ'ffa'Tt^'it'^i.y  6>.!  n- 


3h 


•^Pn(k;t) 


dt 


<   kv  (1  +  ^) 

-   n     v^ 


Since  kv  is  uniformly  boimded  for  all  n  >  1  (see  (91)),  each  term 
of  the  series  (after  the  first  one)  is  of  the  same  order  as  above 
and  the  series  converges  as  well  as  the  one  for  \)/  ,   Consequently 
ij/.  is  continuous  in  (x,y,t)  and  the  second  initial  condition  in  (80) 
is  also  satisfied.   It  might  be  remarked  that  the  estimate  in  (91) 
is  also  sufficient  to  prove  that  \J/  has  continuous  first  partial 
derivatives  with  respect  to  x  and  y,  without  requiring  any  further 
conditions  on  F(k), 

C,    We  now  turn  to  the  stability  question.   Since  the  series  in 
(9i|)  converges  uniformly  with  respect  to  x,y  and  t,  it  is  sufficient 
to  prove  that  each  term  of  the  series  — >  0  for  fixed  y  <  0  and 
X +Ut  as  t  — >  00  ,  Again,  let  N  be  the  lowest  integer  for  which 
[I    "[i  >  0   for  n  >  N,   Then  we  can  write  (8?")  in  the  form 

r  N-1  r       CO   r°    V  ik(v„-U)t 

(95)    ^'  =   EI  J  d^  +  Cl  1   dk^>  A^(k+Ut,y;k)e    ^ 

L  -OO      -^  ^  ^ 

-ik(v^+U)t  1 
+  B^(x+Ut,y;k)e     "     J 

with 

A^(x+Ut,y;k)  =  (ZtiU)  P(k)e        jj-^-p 2v^(^t^-n) 


.-l..,...ik(x.Ut)\(Q>^)^n(^>^)   -n-^ 


„  /  ,  ^  „!.   1  \    CO  tt\-1t,m  ^  ikix+ut;  n  '       n  "  '  n 

B^(x  +  Ut,y,-k)  =  (27IU)  P(k)e         |j^~|j-2 2^^^^^ 
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Consider  first  the  terms  for  n  >  N.   It  has  been  shown  in  Section  I 
that  0  <  (kv^)'  <  U.   Therefore,  neither  (kv^)«  +  U  nor  (kv  )  -  U 
can  vanish  along  the  real   ij^axis.   This  shows  that  k(v  ±  U) 

cannot  have  stationary  points  on  the  path  of  integration.   Then, 

-Ix-e 
since  A  =  0(k  ^   ),  each  term  can  be  integrated  by  parts: 


n 


s 


,00^^  ^ikCv^  -  U)t 


^^^ik(v^  -  U)t^i,  ^  i   ^  -^f-r.:,r-^, .rr^  dk 


"^^  -  -^  J   Kkv^).  -  U]t 

-CD  -co 

and  it  follows  that  the  integral  is  0(t"  ).   The  same  is  true  for 
the  term  containing  B  and  for  all  n  >  N. 

Now  we  will  investigate  the  terms  for  0  <  n  <  N  -  1.   Each  one 
can  be  written  in  the  form  -> 

(96)    "^  c  C  ^2-lZ  eil^f  ^n  -  "'*  .   '"  '  °  e-lXf^n  *   ")*L 


with 


P(k)e^^^^  "^  ^*^  Y^(0;k)  Y^(y;k) 


C^U   +  Ut,y;k)  = -~^- 


n   "^n  "  2 


Since 


2  2-1 
V  ±  U    U   V  g 
n  n^ 


U   -  LL     U  +  V 

V  -  U   ^"  -  ^  v^  +  U 

— has  no  poles  on  the  real   k-axis,  and has  poles  at 

k  =  +  k  ,  where  v  -  U  =  C,   For  the  second  term  in  (96)  x^Je  can, 
therefore,  deform  the  path  L  into  the  real  axis,  and  the  resulting 
integral  is  again  0(t"  )  since  k(v^  +  U)  has  no  stationary  points 
and  we  can  integrate  by  parts. 

To  estimate  the  first  term  in  (96)  we  replace  the  path  L  by 
L  ,  where  L   coincides  with  the  real  axis  except  in  the  neighbor- 
hood of  k  =  +  k  ,  where  it  consists  of  semicircles  of  radius  6  in 
~  n' 

the  lower  half  plane.   The  "phase",  k(v^  -  U),  is  stationary  if 


>   a  >  (. 


.liXW 


(kVjj)«-  U  =  kv'^  -»-(v^  -  U)  =  0 

As  k  -^0,  (kv^).  ->v^(0)    ^7  ^ ^--n 
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-k„  +  k_ 

>- 

L 


and  v^(0)  >  U,   Furthermore,  ■'-'n 

as  |k|  — >  00 ,  (kv  )'  — >  0,   Consequently  there  will  be  at  least 
two  stationary  points,  k  =  ±  k^.   They  cannot  coincide  with 
+  k^,  since  this  would  imply  v^(±kj^)  =  0  which  is  impossible. 
They  cannot  lie  outside  the  interval  (-k^^j  +  k^)  since  v^-  U  <  0 
in  this  region  and  kv^  <  0  for  all  real  k.   Therefore,  they  will 
be  contained  in  ( -k  +5,  k^  -  5 )  for  sufficiently  small  5  >  0, 
The  contribution  to  the  integral  from  this  interval  is  0(t""  Z'^) 
for  some  integer  m  >  2,   The  contribution  from  the  intervals 
mth  |k|  >  k  +  5  is  0(t"  )  as  can  be  seen  upon  integrating  by 
parts,  and  this  leaves  only  the  contributions  along  the  semi- 
circles to  be  accounted  for. 

In  the  neighborhood  of  k  =  +  k^  >  0, 

v^Ck)  =  U  +  (k  -  k^)v^(k^)  +  ... 

For  k'  >  0  we  know  that  V^  <  0;  let  Vj!^(kj^)  =  -t,  o->  0.  Along  the 
semicircle,  k  -  k^  =  5(cos  ©  +  1  sin  9),  ti  <  6  <  2n.  Therefore, 

ik(Vj^  -  U)t    k^5  t(sin  6  -  i  cos  e)  +  0(5^t) 
e  =  e 

and  since  sin  9  <  0  on  the  semicircle,  the  real  part  of  the 

exponent  is  <  0,  and  the  integral  over  the  semicircle  vanishes 

as  t  — ^  00  .   The  same  is  true  for  the  semicircle  near  k  =  -  k^. 

This  completes  the  proof  that  the  solution  ^, given  by  (8?),  tends 

to  a  steady  state  when  observed  from  a  coordinate  system  moving 

mth  the  pressure   disturbance.      It  might   be  noted   that   the 

estimates  for  each  term  in  (8?")  hold  for  general  p.  The  special 

form  for  p  was  used  only  to  show  that  it  suffices  to  consider  each 
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